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. ments on the system size employing a quantum molecular dynamics model. This 

dependence is extracted from the simulations of symmetric collisions of Ca+Ca, 
Ni+Ni, Nb+Nb, Xe+Xe, Er+Er, Au+Au and U+U at incident energies between 
(•^ . 50 A MeV and 1 A GeV. We find that the multiplicity of different fragments scales 

\ with the size of the system which can be parameterized by a simple power law. 
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We report, for the first time, the dependence of the multiplicity of different frag- 



The breaking of colliding nuclei into dozen of fragments has been a central 
topic both in experimental and theoretical Nuclear Physics research. The the- 
oretical interests have been concentrated on the understanding of the reaction 
mechanism behind fragmentation. One has also tempted to look for the role of 
in-medium cross section and nuclear equation of state in fragmentation. These 
studies have been carried out by analyzing the fragment multiplicity and their 
kinetic energy spectra using heavy target/projectile [1-4,6]. However, a very 
little attention has been paid to the system size effects in multi-fragmentation 
although it may not only explore the effect of ratio of surface to volume in 
fragmentation, but can also throw the light on the djTiamical effects which 
are expected to increase with the increasing size of the system. 
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Recently, the FOPI collaboration [3] analyzed the system size effects in frag- 
mentation by studying the reactions of Ni+Ni, Ru+Ru, Xe+Csl and Au+Au. 
They analyzed the spectra by dividing the interacting matter into spectators 
and participants. The spectator fragments show the well known universality [1] 
whereas the intermediate mass fragments emitted from the participant source 
scale with the size of the emitting source indicating the role of the expansion 
of the matter in multi-fragmentation. The role of the expansion has already 
been established in collective flow [7]. 



Motivated by these findings, we present here, for the first time, a complete 
study of the dependence of the fragment multiplicity on the size of the colliding 
(symmetric) nuclei by studying the reactions of different nuclei with masses 
ranging between 40 and 238. We shall show that the multiphcity of all kinds 
of fragments scales with the size of the system which can be parameterized 
by a simple power law oc Al^^{Atot is the mass of the composite system) at all 
incident energies and impact parameters . 



A word of caution should be added here: It has been shown and discussed 
extensively in the literature that the mass yield curve approximately obeys 
a power law behavior oc [2,4]. It has been conjectured (though contro- 
versial) that this behavior is an indication of the phase transition between a 
gaseous and liquid phase of the nuclear matter. This power law behavior ob- 
tained for the mass or charge distribution is for a " given system" . The power 
law dependence, which we are talking about, is something very different. We 
shall show that the multiplicity of a given fragment (or a group of fragments) 
scales with the size of the interacting system which can be parameterized in 
terms of a power law function. 



The basis of our investigation is a quantum molecular dynamical (QMD) model 
[4], where nucleons interact via two and three- bodies which preserves the 
particle correlation and fluctuations. Here, each nucleon is represented by a 
Gaussian wave packet with width y/Z centered around the mean position fi(t) 
and the mean momentum Pi{t): 



The Wigner distribution of a system with {At + Ap) nucleons is given by 
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We shall use here a static soft equation of state along with energy dependent 
nucleon- nucleon cross section. The detail discussion about the different equa- 
tions of state and cross sections can be found in ref.[4]. We follow the nucleons 
till 300 fm/c and then freeze their coordinates. This time is long enough to 
assume that the reaction has finished. The frozen phase space of nucleons 
is, then, clusterized to obtain the fragments. The clusterization is performed 
within Minimum Spanning Tree [MST] method which binds two nucleons in 
a fragment if they are closer than 4 fm. 

Here we simulate several thousand events involving the different symmet- 
ric colliding nuclei like *^Ca Ca,^^ Ni Ni,^^ Nh Nh,^^^ Xe +^^^ 
Xe,^*^^ Er +^'^^ Er,^^'^ Au +^^'^ Au and "^^^U -F^^^ U at incident energies be- 
tween 50A MeV and lA GeV and at different impact parameters h = h/hmax] 
bmax — Ri + R2, Ri is the radius of target/projectile. By varying the size of 
the symmetric system, the system size effect can be studied without varying 
the asymmetry and excitation energy of the reaction. Note that the ALADIN 
experiment varies the asymmetry of a reaction [1] whereas the FOPI experi- 
ments are for symmetric nuclei [3] . The wide range of incident energy between 
50 A MeV and 1 A GeV gives us possibihty to study the multi-fragmentation 
in regions where different reaction mechanisms have been proposed [6]. 

In Fig. 1, we display the largest remnant, the number of emitted nucleons 
and hght fragments with mass = 2 at 50A MeV and 400A MeV. The mul- 
tiplicity of the light mass fragments (LMF's) 2 < A < 4, the medium mass 
fragments (MMF's) 3 < A < 14 and the intermediate mass fragment (IMF's) 
5 <A <mm[Atot/G, 65] is displayed in fig.2. Note that the MMF's and IMF's 
excludes the heaviest fragment. 

The general behavior of all fragments follows the well known trends. In pe- 
ripheral collisions (6 = 0.6), the largest fragment A'^"'^ follows roughly the 
reaction geometry insofar the non overlapping nucleons for the largest rem- 
nant. The size of the heaviest fragment scales with the size of the spectator 
matter which depends on the size of the interacting system. The light charge 
particles ( i.e. the emitted nucleons, the fragment with mass= 2 and LMF's) 
follow the same trend at all incident energies. Their multiplicity is maximum 
at 6 = and which is followed by b = 0.3 and 0.6. The emission of the light 
charge particles depends on the size of the participants which decreases with 
the decrease in the overlapping volume. One also sees that the number of 
light charge particles emitted at 50 A MeV is much smaller than at 400 A 
MeV. Most of the matter at higher incident energies ( i.e. at 1 A GeV) is in 
the form of light charge particles. At higher incident energies, the nucleon- 
nucleon collisions (with large momentum transfer) destroy the correlation in 
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the participant matter and only light particles survive from the reaction zone. 
At lower incident energies, most of the collisions are Pauli blocked and there- 
fore, many nucleons in the reaction zone survive the reaction without suffering 
collisions with a large momentum transfer. The energy received by a target 
in peripheral collisions is not enough to excite the matter far above the Fermi 
level, therefore, a heavier largest fragment A"^"-^, and small number of nucle- 
ons/LMF's and MMF's are emitted at peripheral collisions {b = 0.6). We also 
see that the maximum production of medium mass and intermediate mass 
fragments occurs for 6=0 at 50 A MeV which shifts towards b — 0.3 at 400 A 
MeV [1,4]. With increasing incident energy, the multi- fragmentation becomes 
more and more a phenomenon of peripheral reaction. Due to larger overlap 
between the colliding nuclei at 6 = 0, the excitation in the system at 400 A 
MeV does not allow the formation of heavier fragments, therefore, we see small 
number of MMF's and IMF's. The situation improves aib — 0.3. At peripheral 
collisions {b— 0.6), the amount of excitation, compression and nucleon-nucleon 
collisions decreases and therefore, heavier residual fragment survives. 

The most interesting outcome of figs. 1 and 2 is that in all cases, indepen- 
dent of the mass of the fragment, incident energy ( and excitation energy) 
and impact parameter, the heaviest fragment and the multiplicity of all kind 
of fragments ( i.e. of nucleons, light, medium and intermediate mass frag- 
ments) scales with the size of the system which can be nicely parameterized 
by a power law = c ■ Alg^]Atot is the mass of composite system. The values 
of constants c and r depend on the size of the fragment and on the incident 
energy and impact parameter of the reaction. The dependence of the power 
factor r on incident energy and impact parameter is displayed in figure 3. We 
also tried a function form C-e""^"^, but fits were worse than obtained with a 
power law. It is worth mentioning that a similar power law dependence is also 
obtained for other groups of fragments hke 5 < A < 25. The power law was 
also obtained at other incident energies E= 100, 200, 600 and lOOOA MeV 
indicating the universality of the power law behavior for system size effect in 
multi-fragmentation / clusterization. 

In fig. 3, we plot the value of r as a function of the reduced impact parameter 
b at incident energies 50 and 400 A MeV. The different panels of the figures 
are, respectively, for heaviest fragments free nucleons, fragments with 

mass 2, LMF's, MMF's and IMF's. We see that the value of the r in most of 
the cases depends weakly on the impact parameter. The value of r increases 
with the incident energy and then saturates for very high energies. It is close 
to 2/3 at 50 A MeV whereas it is ~ 1 at 400 A MeV. The only exception in 
this trend is the value of r (= 2.6) for IMF's at 400 A MeV and 6 = 0. This 
exception is rather linked with the mass range used to define an IMF. The 
lower mass limit of an IMF is 5 units and from fig. 1(b), we see that the A^"-^ 
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in this reaction fluctuates around 5 for heavier coUiding nuclei which may lead 
to the emission of few intermediate mass fragments. These fragments cause a 
large value of r. 

We also notice that unlike the disappearance of flow (where the energy of 
vanishing flow varies as A^^^"^) [5], no unique dependence on r exists. This 
behavior of r demonstrates that the degree of clusterization is more for heav- 
ier colliding nuclei which could be attributed to the fast decompression of the 
system in heavier nuclei at central collisions. Similar trends are also reported 
by the FOPI collaboration [3]. If one plots the reported results of FOPI ex- 
periments [3] as a function of size of the system, a similar power- law flt can 
also be obtained [3]. Here one should keep in the mind that the analysis of 
FOPI experiments is done for participant region only. Our present result in- 
cludes both the participant and spectator regions. It is worth mentioning that 
the value of r may depend on the clusterization model one is using. As has 
been noted in Ref.(3), this scaling of the multiplicity of fragments with total 
mass of the system, draws two important aspects: (i) the role of the surface 
to volume ratio in clusterization and second, the dominance of the expansion 
following the compression over the competition between the attractive mean 
fleld and repulsive Coulomb force in multi-fragmentation. The increase in the 
multiplicity of the fragments with mass of the system points toward the fact 
that the expansion of the compressed matter is much faster in heavier systems 
compared to lighter systems. 

Summarizing, using quantum molecular dynamics (QMD) model, we flnd, for 
the flrst time, a power- law behavior for the system- size effect in fragmen- 
tation at all incident energies between 50 A MeV and 1 A GeV and for all 
colliding geometries. This dependence exists for any kind of fragment/group 
of fragments and can be parameterized in terms of a power law. 
This work was supported by Council of Scientific and Industrial Research grant 
No. 03 (0823) /97/EMR-II. 
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Fig. 1. The heaviest fragment A'^'^^ , the multiphcity of free- nucleons and fragments 
with mass A = 2 as a function of composite mass of the system Atoti= + ^p) 
at different scaled impact parameters h. The left hand side is at 50A MeV whereas 
right hand side is at 400A MeV. The displayed results are at 300 fm/c. 
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Fig. 2; " Study of system- size elfecis " by J. Singh et al. 



Fig. 2. Same as fig. 1, but for light mass fragments (LMF's) {2 < A < 4), medium 
mass fragments (MMF's) (3 < ^ < 14) and intermediate mass fragments (IMF's) 
(5 < A < mm[l/6^toi,65]). 
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Fig. 3: " Study of system- size effects " by J. Singh et al. 



Fig. 3. The r as a function of the reduced impact parameter h. The solid and dashed 
lines shows, respectively, the results at 50 A MeV and 400 A MeV. 
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